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I. INTRODUCTION 


CP violation [CPV) was first seen in the decay of neutral kaons Q. It was shown some time ago that CPV is a 
necessary but insufficient constraint in order to generate a net baryon anti-baryon asymmetry in the universe [^. In 
the ensuing years there has been a tremendous amount of activity by the high energy physics community to better 
understand the role of CPV in our model of nature — the Standard Model of Particle Physics (SM). 

CPV in the SM is described by a single complex phase in a 3 x 3 quark-mixing matrix, Vckm, called the CKM Q 
matrix: 
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Vckm describes the coupling of the u, c and t quarks to d, s and b quarks, which is mediated by the exchange of a 
W boson. In B-meson decays the interesting CP violating parameters of the SM are related to the angles (/3, a, and 
7 )^ and sides of the so-called Unitarity Triangle (as shown in Figure QJ. The angles are defined as: 


a = arg[-y,U*,/P,,U„*,], 

(2) 

/3 = arg[-U,dU*,/y,y;], 

(3) 

^ = ai:g[-V^^V:jVM]- 

(4) 



FIG. 1: The Unitarity Triangle. 


This review is a summary of the experimental constraints on the Unitarity Triangle angle a obtained from 
decays involving b ^ u transitions. The CKM angle f3 is well known and consistent with SM predictions 
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^ The Belle Collaboration use the notation (pi, <p 2 , and cps for the angles of the Unitarity Triangle. 
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constraint on a constitutes a test of the SM description of quark mixing and CPV in B-meson decays. A significant 
deviation from SM expectation would be a clear indication of new physics. It is possible to obtain a SM prediction 
of a from indirect constraints by combining measurements of the CKM matrix elements |t4is|, |Kid|, IKt&l, and \Vcb\^ 
CPV in mixing from neutral kaons, B-B mixing in Bd and Bs mesons and the measurement of sin 2/3 from b —> ccs 
decays. The SM predictions for a are (98.2 ± 7.7)° from and ( 97 ^ 19 )° @- 

Measurements of a have recently been performed at the B-factories, the BABAR experiment 0 at SLAG and the Belle 
experiment Q at KEK. The B-factories study the decays of B-mesons produced in the process e+e“ ^ T(4S) ^ BB, 
where BB is either B+B“ or B°B . The neutral B-mesons are produced in a correlated P-wave state. Until one 
of the B-mesons in an event decays, there is exactly one B° and one B^ meson. The main goals of the B-factories 
include the study of CPV in B-meson decay, and to over-constrain the position of the apex of the Unitarity Triangle 
{p,r]). The final states of interest for the study of a are B ^ hh' , where h = p,7r. Figure |21 shows the Feynman 
diagrams corresponding to the dominant amplitudes contributing to the h'^h'~ decays. 



FIG. 2: The tree (left) and gluonic loop (right) contributions to B ^ h' decays. 


Interference between the amplitude for direct decay, and decay after B°-B° mixing to a h^h'~ final state, results in 
a time-dependent decay-rate asymmetry between B° and B° decays that is sensitive to the CKM angle a. To study 
the time-dependent asymmetry with neutral B-mesons, one needs to measure the proper time difference At between 
the decay of the two B-mesons in the event. On reconstructing the h'^h'~ final state of interest (Bjec), one needs to 
determine the flavor of the other B-meson (Btag)^ from its decay products. The time difference between the decays of 
the two neutral B-mesons in the event (Brec, Bt^g) is calculated from the measured separation Az between the B^ec 
and Btag decay vertices. The B^ec decay vertex is determined from the two charged-pion tracks in its final state. The 
Btag decay vertex is obtained by fitting the other tracks in the event. 

The signal decay-rate distribution of a CP-eigenstate decay, /+(/_) for Btag= B°(B°), is given by: 

g-AdZ-r 

f±{At) = —--[1 ± S sin(AmdAt) =F C cos(AmdAt)], (5) 

4t 

where r = 1.536 ± 0.014ps is the mean B° lifetime and Arrid = 0.502 ± 0.007ps“^ is the B°-B° mixing frequency j^. 
This assumes that there is no difference between B° lifetimes, AP = 0. The parameters S and C are defined as: 


2/mA 

^ ^ 1 + |A|2’ 


^ 1-|AP 
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( 6 ) 


where A = is related to the level of B*^-B° mixing (q/p), and the ratio of amplitudes of the decay of a B° or B° 
to the final state under study {A/A). CPV is probed by studying the time-dependent decay-rate asymmetry 


R{At)-R{At) 

R{At)+R{At)' 

where R{R) is the decay-rate for B°(B°) tagged events. This asymmetry has the form 

A = S sm{AmdAt) — Ccos{AmdAt). 

Belle use a different convention to BABAR with C = —Acp- 


( 7 ) 

( 8 ) 


^ i.e. determine if Stag is a or . 
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In this article, all results are quoted using the S and C convention. In the case of charged S-meson decays (and 
7 r*^ 7 r*^ as there is no vertex information) one can study a time integrated charge asymmetry 


. N-N 

= WTJi’ 


(9) 


where N {N) is the number of B (B) decays to the final state. A non-zero measurement of S, C or Acp for any of 
the decays under study would be a clear indication of CPV. 

In the absence of loop (penguin) contributions in decays to the CP eigenstate h~^h~, S = sin 2a, and (7 = 0. 
The presence of penguin contributions with different weak phases to the tree level amplitude shift the experimentally 
measurable parameter agg away from the value of a. In the presence of penguin contributions a^s = a + (^ajpgnguin), 
S = sin2aeff, and C can be non zero. 

For decays to p^p~, S= Sp or St and C= Cl or Ct are the CP asymmetry parameters for the longitudinal 
and transversely polarized signal, respectively. The decay-rate distribution for decays to the p^ir^ final state is 
more complicated than the above description as is not a CP eigenstate (discussed below). 

There are two classes of measurement that the B-factories are pursuing. The main goal is to use SU{2) isospin 
relations to relate different hh' final states and limit 5a(penguin) in each of the modes ^ tt+tt”, 5° ^ and 

B^ —> p'^p~. In some cases only weak constraints on a are obtained when using SU(2) relations and one must use a 
model dependent approach to obtain a significant result. 

The remainder of this article describes the isospin analysis used to constrain (5a(penguin)! and experimental tech¬ 
niques. There is a discussion of the results at the end of the article. Throughout this article, the experimental results 
are quoted with statistical errors preceding systematic errors unless otherwise stated. 


II. ISOSPIN ANALYSIS OF B ^ hh' 


One can use SU{2) isospin to relate the amplitudes of B decays to tttt final states This results in two relations: 


-^A+~ = A+° - A°°, 

V2 

J-A^- = A~° 

V2 


( 10 ) 

( 11 ) 


where A^ are the amplitudes of B (B) decays to the final state with charge ij. These two relations correspond 
to triangles in a complex plane as shown in Figured 



FIG. 3: The isospin triangle for B — > tttt decays. 

This approach only considers tree and gluonic penguin contributions. Possible contributions from electroweak 
penguins (EWP) are ignored as they do not obey SU{2) isospin symmetry. EWPs have the same topology as the 
gluonic penguin diagram in Figure|21 with the gluon replaced by 7 or bosons. In the absence of EWP contributions 
|A+°| = i.e. Acp = 0 for B+ —> 7 r+ 7 r°. After aligning and the phase difference between A+ and 

A is 2(5a(penguin)- Order to measure a one must measure the branching fractions (B) and charge asymmetries 
(Acp, or C) of B decays to tt+tt", 7 r° 7 r°. 
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The decays B pp are more complicated; one has an isospin triangle relation for each of the three amplitudes in 
the hnal state, one corresponding to the CP even longitudinal polarization and two corresponding to the CP even, 
and CP odd parts of the transverse polarization. The p^ir^ final state is not a CP eigenstate. This results in the 
need for a pentagon isospin analysis ^ Dalitz Plot (DP) analysis of the 7r+7r“7r° final state 0 to measure a. 

There are several assumptions implicitly used in the isospin-based direct measurements of a. The assumptions 
common to all decay modes are (i) to neglect EWP contributions and (ii) to neglect other SU{2) symmetry breaking 
effects. In addition to this, possible 1=1 amplitudes ^ in P^P~ decays are ignored. Several mou ps have 

estimated the correction due to the SU(2) breaking effect of EWP contributions to be 0(1.5 — 2)° [ij,[l3- These 
estimates consider contributions from the two EWP operators assumed to be dominant in the effective Hamiltonian. 
Some estimates of other expected SU{2) symmetry breaking effects exist, where any correction is estimated to be 
much less than the current experimental precision [l^ . 

III. EXPERIMENTAL TECHNIQUES 

Continuum e+e” ^ qq (q = u, d, s, c) events are the dominant background to B ^ hh' decays. Signal 
candidates are identihed using two kinematic variables: the difference between the energy of the B candidate and 
the beam energy ^/s/2 in the center of mass (CM) frame AE; and the beam-energy substituted mass mss = 
yj (s/2 -I- Pi • TpsYlE'f — p^. As the H-factories collect data at the r(4S) resonance y/s = 10.58 GeV. The B 
momentum ps and four-momentum of the initial state (EijPi) are defined in the laboratory frame. Event shape 
variables are used to further discriminate between signal and continuum background |17| . 

The flavor of Stag is determined from its final state particles. Both BABAR and Belle characterize the final state 
particles of i?tag into events where there are leptons, kaons and 7r-mesons. The charge of these reconstructed tracks 
in the hnal state indicate if Stag is a B^ or a 5°. The signal purity varies depending on the momentum and type of 
tracks in the hnal state. Events with leptons in the hnal state are the cleanest and those with charged kaons in the 
hnal state have less background than those with 7r-mesons. In addition to this, the probability of assigning the wrong 
havor to Htag increases as the background increases. BABAR separates events into mutually exclusive categories of 
events, whereas the Belle combines this information into a single variable r, and tunes event selection in bins of r. 
The two approaches are equivalent and are described in more detail in the following references [ 30 . 

The At distributions of Equations 13 and M are convoluted with a detector resolution description, which differs for 
signal and continuum background jl fll l2(l | , and also take into account dilution from incorrectly assigning the havor 
of Htag. The signal and background parameters are extracted using an extended unbinned maximum-likelihood ht to 
the data for the analyses described here. 


IV. RESULTS 

A. B ^ TTTT 

The simplest de cay s to study in the pursuit of a are B —> tttt. Both experiments have measured B —>■ 7r'''7r“, 
and 7r°7r° decays |^. The results are summarized in the Table below. The BABAR measurements use a data 
sample of 210/6“^, the Belle result use 78fb~^, and the other Belle results use 253/6“^ of integrated luminosity. 
All of these modes are now well established experimentally and so it is possible to perform an isospin analysis. The 
charge asymmetry measurement of B^ -i- 7r°7r° accounts for the effect of B^-B^ mixing. 


mode 

Expt. e(xio-®) 

S 

C 

7r'’'7r~ 

BABAR 5.5 ± 0.4 ± 0.3 

Belle 4.4 ±0.6 ±0.3 

-0.30 ±0.17 ±0.03 

-0.67 ±0.16 ±0.06 

-0.09 ±0.15 ±0.04 

-0.56 ±0.12 ±0.06 


mode 

Expt. 

B (xlO-®) 

Acp 

7r^7r° 

BABAR 

5.8 ±0.6 ±0.4 

-0.01 ±0.10 ±0.02 


Belle 

5.0 ±1.2 ±0.5 

-0.02 ±0.10 ±0.01 

o 

o 

BABAR 1.17 ±0.32 ±0.10 

-0.12 ±0.56 ±0.06 


Belle 

9 q+ 0.4 +0.2 
^•o_o.5 _o.3 

-0.44t° g ± 0.17 


One should note that the Belle measurement of B^ —> tt^tt constitutes an observation of CPV in this decay at 
a level of 5.4cr, and evidence for direct CPV at a level of 4 .Oct. This is the second observation of CPV in B-meson 
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decays, the first being the measurement of a non-zero value for sin 2/3 from b ccs decays 0. The BABAR data 
are consistent with the hypothesis of no CPV in this decay. In recent years the two sets of results have started to 
converge to a common value, but more statistics are required to resolve the controversy. 

The TTTT isospin analysis is limited by the value of B{B^ —> 7r°7r°). Unfortunately this branching fraction is nei¬ 
ther large enough to provide sufficient statistics with the current data set to perform a precision measurement of 
(5a(penguin)j HOI' Small enough to enable a strong bound on this quantity. BABAR have performed an isospin analysis 
resulting in |^Q;(penguin)| < 35° (90% C.L.). Belle’s data exclude values of a between 19 and 72° (95.4% C.L.), and 
constrain l^a/penguin)! < 38° (95.4% C.L.) |^. One requires a significant increase in statistics to perform a precision 
measurement of a using B tttt decays. However, it is possible to extract a model dependent constraint on a from 
the current results assuming SU(3) symmetry (i.e. exchanging u and s quarks), Gronau and Rosner obtain a value 
of a = (107± 13)° m. 


B. B ^ pp 

The decay B pp is that of a spin zero particle decaying into two spin one particles (as shown in Figure As a 
result, the CP analysis of B decays to p^p~ is complicated by the presence of three helicity states {H = 0,±1). The 
H — 0 state corresponds to longitudinal polarization and is CP-even, while neither the H = +1 nor the H = —1 state 
is an eigenstate of CP. The longitudinal polarization fraction is defined as the fraction of the helicity zero state in 
the decay. The angular distribution is 

^ (II cos^ 01 cos^ 02 + ^(1 - /l) sin^ 0i sin^ 92 ] (12) 

Fd cos a cos ^2 4 \ 4 J 

where 0i (i = 1, 2) is defined for each p-meson as the angle between the tt^ momentum in the p rest frame and the 
flight direction of the P° in this frame. The angle 4> between the p-decay planes is integrated over to simplify the 
analysis. A full angular analysis of the decays is needed in order to separate the definite CP contributions of the 
transverse polarization; if however a single CP channel dominates the decay (which has been experimentally verihed), 
this is not necessary |^ . 


3 



FIG. 4: A schematic diagram of the decay of a P-meson via two vector particles, Vi and V 2 , into a four-particle final state. 

The longitudinal polarization dominates this decay 111123 . Not all of the pp final states have been observed, 
however as B{B^ P^) is small, one can conservatively assume that it is longitudinally polarized when performing 
an isospin analysis. The measured branching fractions, /i, and C for P —> pp are summarized in the Table below, 
and Sl = -0.33±0.24+°;)'® (0.09±0.42±0.08) for BABAR. (Belle). The BABAR (Belle) p+p" results use 210 (253)/5-i 
of integrated luminosity, respectively^. The BABAR (Belle) p'^p^ results use 82 (78) fb~^, and the upper limit on p°p° 


^ The branching fraction reported by BABAR is from 82 fb 
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mode 

Expt. 

B (xlQ-®) 

U 

Cl 


BABAR 

Belle 

30 ± 5 ± 4 
24.4 ± 2.2t® ® 

0.978 ± 0.014t° °2i 

n +0-033 +0.029 

U.yOi_o.o39 -0.031 

-0.03 ±0.18 ±0.09 
0.00 ± 0.30tg;J^ 


mode 

Expt. 

B (xlO-®) 

/l Acp 


BABAR 

22.5tl:l ± 5.8 

0.97tg;2? ± 0.04 -0.19 ± 0.23 ± 0.03 


Belle 

31.7±7.lt^;f 

0.95 ±0.11 ±0.02 0.00 ±0.22 ±0.03 


BABAR 

< 1.1 (90% C.L.) 

- 


uses 210/6“^ of integrated luminosity. Recent measurements of the 5+ ^ branching fraction and u ppe r limit 

for 13 indicate small penguin contributions in R ^ pp, as predicted by some calculations |3 • The 

ultimate uncertainty on a from B —> pp will depend on the branching fraction and CP content of —> p^p^. The 
decay B^ p^p^ has an all charged track final state and it is possible to measure both S and C for the different 
CP-eigenstate components of the decay, unlike B^ 7r°7r° decays where one can only measure C. 

Given that penguin pollution is small, it is possible to perform an isospin analysis of the longitudinal polarization 
of the B ^ pp decays, and use the results of BABAR's time-dependent CP analysis of p^p~ |3 to constrain a. If 
one does this, using the aforementioned assumptions, one obtains a = (100 ± 13)° using BABAR data |3- The error 
on a is dominated by |<5a(penguin)| < 11° (68% C.L.). Belle have recently produced a preliminary measurement of 
the polarization, branching fraction and CP parameters of ^ P~^P~ The Belle constraint is a = (87 ± 17)°, 
which is slightly weaker than the BABAR result. Both pp isospin analyses use the same experimental information for 
the decays ^ and B^ FigureElshows the (1-C.L.) plot of a corresponding to the isospin analysis 

of the longitudinally polarized pp data [^. There is also a non-SM solution for a near 175°. 

Our knowledge of ^Q;(penguin) is primarily determined from the experimental knowledge of the p^p^ and p^p^ 
branching ratios. As in the case of P —*■ tttt, our knowledge on a is limited by the penguin contribution. However, 
for B pp decays (5a(penguin) is sufficiently small to allow the P-factories to perform a meaningful measurement. 



a (deg) 


FIG. 5: A plot of (1 - C.L.) of a determined from the BABAR (dashed) and Belle (dotted) pp isospin analyses 0- The shaded 
region corresponds to the constraint obtained from the average of the two experiments. The point with error bars corresponds 
to the la interval (68 % C.L.) obtained from the CKM fit when excluding direct measurements of a. 

The pp isospin triangle does not close when taking into account the measurements of the P-factories. Both the 
CKM Fitter and UTfit groups have noted that SU{2) symmetry constraints prefer a smaller branching fraction for 
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the decay, at the level of 1.5 — 2(t. The measurements did not test for possible non-resonant backgrounds, 

S-wave TTTT contributions under the or account for exclusive B backgrounds to decays excluding charm particles. 
It is important to see that this decay mode is studied in more detail with higher statistics in the future. 


C. B ^ p-K 


The decays B ^ pir can be analysed in two different ways. The more straightforward approach is to cut away 
interference regions of the 7 r+ 7 r“ 7 r° DP and analyze the regions in the vicinity of the p resonances. This is the so-called 
Quasi-2-body approach (Q2B) and it avoids the need to understand the interference regions. The drawbacks of the 
Q2B method are that one looses a lot of statistical power by cutting on the DP, and that interference is unaccounted 
for. A corollary of this is that one requires more statistics than the ^-factories currently have in order to obtain a 
significant constraint from the pentagon isospin analysis of i? —> p 7 r. The alternative is to perform an analysis of 
the B —> TT+TT^TT*^ DP, accounting for the interference between intersecti^ presonance bands and other resonant 
structures. The Q2B approach has been studied by BABAR and Belle |^, and the first time-dependent DP 
analysis has been performed by BABAR [^. 

In the Q2B approach, one fits a time-dependence of 

-\At\/T 

f±{At) = (1 ± Acp )— - - [{S ± AS sin(A 77 i(iAt) — (C ± AC) cos(AmrfAt)], 

4t 

(13) 

where there are three additional parameters in comparison with Equation 0 These parameters are a charged asym¬ 
metry, Acp, between decays to p+ 7 r“ and p“ 7 r+ final states, and two dilution parameters; AS and AC. The details of 
the DP analysis can be found in |34| . where one varies a larger number of parameters in the nominal fit, and converts 
these to the same observables as the Q2B approach. The Table below summarizes the experimental constraints on 
B ^ pTT (DP analysis from BABAR and Q2B analysis from Belle). The branching fractions measured for this decay 
are (22.6 ± 1.8 ± 2.2) x 10“® and ( 20 . 8 ^ 6,3 I 3 . 1 ) x 10“® by BABAR and Belle, respectively [s^. These results used 82 
(29.4)/5“^ of integrated luminosity and the time dependent results use 192 (140)/6“^, respectively. 


Expt. 

BABAR 

Belle 


S C Acp 

-0.10 ± 0.14 ± 0.04 0.34 ± 0.11 ± 0.05 -0.088 ± 0.049 ± 0.013 
-0.28 ± 0.23lo;o8 0.25 ±0.1712 Kg -0.16 ± 0.10 ± 0.02 


Using BABAR’s DP result, one obtains the following constraint; a = (113l37±6)°. FigureElshows the corresponding 
(1 - C.L.) plot for a. This result is self-consistent as the strong phase differences and amplitudes are determined solely 
from the structure of the DP. The unique aspect of this result is that there is only a single solution between 0 and 
180°, and thus a two fold ambiguity on a. As a result this measurement is an important constraint on a, ancillary to 
that from B ^ pp. 

It is possible to derive a model dependent constraint on a from these data as shown by Gronau and Zupan |36| . 
One can make several assumptions, including (i) the relative strong phase differences between tree level contributions 
for p+ 7 r“ and p“ 7 r+ being less than 90° (ii) SU{3) being exact for penguin amplitudes (hi) neglecting EWP (iv) 
neglecting annihilation and exchange diagrams to derive a constraint on a. SU{3) breaking effects are estimated in 
this model using CLEO, BABAR and Belle data. When doing this Belle obtains a = (I02± llg^pti 15moiiei)°- Gronau 
and Zupan compute the corresponding result for the BABAR analysis as (93 ± ± 15modej)°- There are other 

models proposed in the literature, for example 


V. DISCUSSION 

The indirect measurements of a from the SM are (98.2 ± 7.7)° i and (97^21)° @ from the UTfit and CKM Fitter 
groups, respectively. Direct measurements from the B-factories are in good agreement with these predictions, where 
the uncertainty on a from B —> hh' decays is known to 0(9°). The precision of this result is dominated by an 
isospin analysis of B —> pp decays, with an important contribution from B —> p 7 r as the latter result suppresses a 
second solution for a near 175°. This solution is disfavored at the 1.3cr level by studies of B° ^ 7 r+ 7 r“ 7 r° decays. 
The B-factory constraints on a will be improved as additional data are analyzed over the next few years, and these 
measurements will start to test the closure of Unitarity Triangle with precision. A significant isospin analysis result 






a (deg) 


FIG. 6: A plot of (1 - C.L.) of a determined from BABAR’s tt® DP analysis 1^ . 


from B TTTT requires a considerable increase in statistics. The measurement of a from B —>■ tttt remains a vital part 
of the program, as theoretical uncertainties in the extraction of a (other than the determination of ^Q;(penguin)) are 
the best understood of all the channels that have been discussed here. 

Similarly, the model dependent calculations on a are in agreement with the SM. The constraints from B — > tttt 
and p'K decays have a precision of 13° and 15.5°, respectively. These constraints are dominated by theoretical 
uncertainties. The precision on the model dependent constraints obtained are comparable with the isospin analysis 
result from B pp. 

The measurement of a using B —> oitt decays was proposed some time ago j^. The recent observation of 
B^ a^7r“ with a branching fraction 0(45 x 10“®) raises the question of how well one can ultimately measure 
a using this decay mode. The use of SU{3) to measure a with B oitt decays has recently been proposed [iof. 

Studies of the LHCb experiment’s sensitivity to a using B 7r+7r“7r® decays have been performed LHCb 
is expected to achieve a 10° precision on a within the first year of data taking. The physics reach of decays involving 
b ^ u transitions with all charged particles in the final states, such as B^ and B^ a^7r“ should also 

be investigated. The current i?-factories will reach a precision of 7 — 10° on a with 1 ab~^ of data using B ^ pp 
decays. The ultimate precision reached depends on the values of the branching fraction, Jl, S and C in ^ 
Sensitivity projections for the 7r'''7r“7r® DP analysis are very dependent on the many input parameters involved. For 
this reason, the running experiments have not extrapolated the precision of a obtained from 7r“''7r“7r® decays to higher 
luminosities. 
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